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Abstract 

Using a coherence time extracted from high precision proton-nucleus Drell-Yan mea- 
surements and a nuclear absorption cross section extracted from pA charmonium produc- 
tion experiments, we study J /tp production and absorption in nucleus-nucleus collisions. 
We find that coherence time effects are large enough to affect the measured J/^-to-Drell- 
Yan ratio. The S+U data at 200A GeV/c measured by NA38 are reproduced quantita- 
tively without the introduction of any new parameters. However, when compared with 
recent NA50 measurements for Pb+Pb at 158A GeV/c, the data is not reproduced in 
trend or in magnitude. 
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Ultrarelativistic heavy ion collisions offer the tantalizing possibility of forming and studying 
a new form of matter predicted by QCD: the quark-gluon plasma. A vigorous experimental 
program has existed at the CERN SPS for more than ten years. RHIC now signals the dawn 
of a new era in heavy ion physics at Brookhaven National Laboratory. Several experimental 
signals have been put forward as candidates for QCD plasma signatures []!]]. Of those, the most 
famous is probably that of J /if) suppression in nucleus-nucleus collisions. The theoretical and 
experimental activity that have followed this seminal suggestion have been considerable as the 
disappearance of the J /if) can directly be linked to deconfinement and Debye screening in the 
plasma ||. The interested reader will find a recent snapshot of the state of this field in Ref. 

I- 

Before an experimentally observed J / if) suppression pattern is interpreted as an unambigu- 
ous signal of the existence of a quark-gluon plasma, it is imperative to rule out all competing 
explanations of purely hadronic origin. Moreover, the hadronic scenarios considered should 
incorporate elements of physics that are known to be relevant at the energy scale under con- 
sideration. It is one such line of thought that we follow in this paper. We study charmonium 
production in relativistic heavy ion collisions along with the appropriate background, Drell- 
Yan pair production. Our paper is organized as follows: First we recall the main features of a 
model that is successful in explaining high precision Drell-Yan data measured in proton-nucleus 
collisions. Those data enable one to extract a formation time characteristic of the emission of 
soft hadrons, essentially pions. Next we recall the application of this model to the produc- 
tion of J / if) in pA collisions. From those measurements we have extracted a cross section for 
J /if) absorption on the nucleon. With this formulation we make parameter-free calculations for 
nucleus-nucleus collisions and compare them with experimental data. 

In almost all considerations involving heavy ion collisions at any energy, the issues of dy- 
namics and elementary processes remain intimately connected and inseparable. In view of this, 
a successful modeling of nuclear collisions is a necessary prerequisite to a deeper exploration of 
the fine points of the nuclear dynamics. To simulate the heavy ion collision we prefer to work 
with hadronic variables rather than partonic ones, and make a straightforward linear extrapo- 
lation from proton-proton scattering. This extrapolation, referred to as LEXUS, was detailed 
and applied to nucleus-nucleus collisions at beam energies of several hundred GeV per nucleon 
in Ref. Q. Briefly, the inclusive distribution in rapidity y of the beam proton in an elementary 
proton-nucleon collision is parameterized rather well by 

W l (y) = + (1 - X)S(y - y) , (1) 

smh y 

where yo is the beam rapidity in the lab frame. The parameter A has the value 0.6 independent 
of beam energy, at least in the range in which it has been measured, which is 12 — 400 GeV ||. 
It may be interpreted as the fraction of all collisions which are neither diffractive nor elastic. 
As a nucleon cascades through the nucleus its energy is degraded. An underlying assumption 
in this model is that of straight line trajectories. In the case of a nucleus-nucleus collision, one 
obtains the single-particle rapidity distribution of the m'th projectile nucleon after a collision 
with the n'th target nucleon through the solution of an evolution equation M: 



(y) = J ^vp^vt^^u-MpW^-i^vq - vt) Q{y - vt, vp -yr,y- yp) ■ ( 2 ) 
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In the above, the kernel is 



posh 

Q(s,t,u) = X^-^ + (l-X)5(u), (3) 

originating from Eq. (|l|). Equation (|]) is a Boltzmann-like equation that is solved numerically. 
This rapidity distribution then gets folded with impact parameter over the density distributions 
of the projectile and target nuclei, using a method described in detail in Ref. ||. 

Recently we have extracted the quantum coherence time needed to reproduce Drell-Yan 
pair production data in pA collisions. This can also be formulated in terms of the Landau- 
Pomeranchuk-Migdal effect ||. We briefly recall the procedure followed. We began by com- 
puting the Drell-Yan yield at leading-order (LO) with the GRV structure functions f?j with a 
K factor. Those structure functions reflect a flavor-asymmetric Dirac sea. Adopting a fixed K 
factor of 2.1, we compared the results to pp Drell-Yan data at 800 GeV/c || and found the 
agreement to be excellent for all measured values of Xp. 

Turning then to the case of proton-nucleus collisions as measured by the E772 collaboration 
0, we deduced [JTOj that the formation (or coherence) time needed to fit the measured a^X /Vpo 
ratios at different values of xp is r = 0.4 ± 0.1 fm/c, in the frame of the colliding nucleons. 
Practically, this coherence time can be related to an initial state energy loss for some of the 



Drell-Yan producing collisions [TT| as follows. In LEXUS, we assumed that the energy available 
to produce a Drell-Yan pair was that which the proton had after n previous collisions. In order 
to reproduce the 800 GeV/c E772 data, we needed n = 5 ± 1. The n collisions correspond to 
a path length of n/ajffxp in the target nucleus rest frame, where a^ N is a total cross section 
of 40 mb, and p is a nuclear matter density of 0.155 nucleons/fm 3 . Lorentz-transforming to 
the nucleon-nucleon center of mass, one obtains the value of the proper coherence time quoted 
above. In this language, a traditional Glauber-type model (with no energy loss) would have 
n = oo. Fixing r , the range in n appropriate for the energies being considered in this work 
(nucleon momenta of 158 GeV/c and 200 GeV/c) is found to be 2 < n < 3. 

We then investigated J/tp production in pA collisions [|T2|, |T3J. The additional input needed 
there was the cross section for producing J / ip in elementary nucleon-nucleon interactions. We 
used a parametrization that follows from a tabulation of data by Lourengo ]T3| : 

Ba NN ^j/^(x F > 0) = 37(1 - mj/^/i/s) 12 nb . (4) 

Here, yfs in the center of mass energy of the nucleon pair and B is the branching ratio into a 
muon pair. Using the functional xp dependence of the differential cross section as measured by 
E789 [fbl, one can write a normalized differential cross section to use as an input in LEXUS: 



dGN ^ = Qa NN ^j/4xF > 0)(1 - M) 5 . (5) 
axp 



From our analysis |12[], we have extracted a J/ip absorption cross section in nuclear matter 



of 3.6 mb. It is worthwhile to note that this value is in numerical agreement with the same 
quantity deduced from experiments of J/tp photoproduction on nuclei [15]. It is smaller than 



that used in other phenomenological heavy ion applications jllj. The parameters in this model 
are thus completely determined by proton-nucleus data. 



We now turn to recent experiments on the production of the J/ip in S+U [|17| collisions 



at 200A GeV/c and in Pb+Pb fL8| collisions at 158A GeV/c, at the CERN SPS. Because the 
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Figure 1: Ratio of J/ip to Drell-Yan production cross sections. B is the branching ratio into a 
muon pair. The data is from Ref. |17] . 



J/tp is measured through its decay into dimuons, the production cross section has traditionally 
been divided by the natural background in the appropriate invariant mass region: that of Drell- 
Yan pairs. However, since the absolute cross section measurements are now available, we will 
first verify the predictions of our model there. Including the appropriate respective detector 
acceptance, the results for absolute cross sections are show in Table [[]. 

Consider the system S+U at 200A GeV/c. The coherence time arguments made earlier 
in this paper suggest that the values n = 2 and n = 3 should bracket the experimental data. 
We observe that indeed this is so, both for the measured Drell-Yan and J/ip absolute cross 
sections. One can go further and plot Ba J ^ / a DY against collision centrality. One needs a 
model to map the impact parameter bins that enter as input in our dynamical model into bins 
of measured transverse energy. The experimental collaboration has in fact provided the impact 
parameter range that corresponds to a measured J / ^-Drell-Yan ratio ]17| . Comparison to the 
experimental data is shown in Fig. |I[ One can see that our results are consistent with the 
data within experimental uncertainties. Again, we emphasize that no new parameters were 
introduced. It is also worthwhile to note that the numerator and denominator of the plotted 
ratio have been calculated from "first principles" , the meaning of which is clear in the context of 
this paper: LO Drell-Yan calculations and a parametrization of the differential J/ip production 
in nucleon-nucleon collisions. The absolute cross sections calculated with no energy loss (or 
infinite coherence time) fail to reproduce the experimental results. This is also the case for 
their ratio. 

We now turn to experimental results obtained by the NA50 collaboration with Pb projectiles 
and targets. From Table p] we see that the measured Drell-Yan cross section exceeds our larger 
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Figure 2: Same as Fig. [I], except that this is for the heavier system. The data is from Ref. ||19|| . 



(n = 3) value by about one standard deviation. The experimental J/tp value falls within the 
predicted range. Plotting the J/^-to-Drell-Yan ratio against the impact parameters determined 
by the experimental collaboration |1| one obtains Fig. |2|. Application of this model with its 
parameters determined solely from pA physics does not yield a satisfactory representation of 
this experimental data. The latter is not reproduced in trend nor is it in magnitude. Note, 
however, that the poor quality of this fit is entirely comparable with those obtained with other 
hadronic approaches p0 |. Also shown in this figure is the effect of the coherence time on 
this ratio. It appears that this effect is not as spectacular here as it partially cancels in the 
numerator and denominator. The flattening and slight increase of the ratio, as one goes to 
smaller impact parameters, can be attributed to the J /if) cross section growing slightly faster 
than the Drell-Yan. Here also, the calculated absolute cross sections with no energy loss far 
exceed the experimental values. 

We also considered possible nuclear structure effects on the ratio shown in Fig. [2| It is 
known that parton distribution functions that have a flavor-asymmetric Dirac sea, like the one 
we use in this work, will yield different Drell-Yan cross sections depending on whether one has 
p + p, p + n, n + norn + p collisions. In our treatment the isospin content of the nucleus is 
assumed to be uniformly distributed according to the overall charge of the colliding partners. 



Experimentally, Pb is known to have a neutron skin |2l] of 0.19 ± 0.09 fm. This value is in 
fact too small to have an effect on the calculations shown here. Finally, it seems useful to point 
out that in Fig. the Pb data does not seem to converge to the vacuum ratio as one moves 
towards more peripheral collisions, unlike the measurements of S-induced reactions. 
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We have investigated nucleus-nucleus collisions with a model that incorporates the coher- 
ence time associated with the emission of soft quanta in hadronic interactions. This approach 
translates into lost energy for the formation of hard radiation, such as high-mass Drell-Yan 
pairs and J/if>. We have obtained results in quantitative agreement with experimental data for 
the reaction S on U at 200A GeV/c. The ratio of J/ip to Drell-Yan cross sections as a function 
of collision centrality, as well as the total absolute cross sections are reproduced by our model. 
Therefore, we can understand Drell-Yan and J/ip formation in pA and S+U collisions in terms 
of the same physics. This model fails to reproduce measurements done in connection with the 
heavier Pb+Pb system. 

Several points still need to be clarified. It will be very instructive to repeat this analysis in 



partonic variables including nuclear shadowing [22]. A systematic exploration of the freedom 
allowed by the most recent high-precision pA measurements |S| is called for and is underway. 
Nevertheless, if the Pb+Pb data stand the test of time, it does not seem possible to escape 
the conclusion that J/ip suppression is caused by high energy density. Whether it is due to 



absorption on hadronic co-movers |24] or quark-gluon plasma remains an open and exciting 
question. 
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Table 1: The S+U processes are measured by NA38 |17| in collisions at 200A GeV/c and 
the Pb+Pb processes are measured by NA50 [|TB[ in collisions at 158A GeV/c. The entries in 
columns labeled with different values of n represent calculated values. The J j ip cross sections 
are obtained using a J/^-nucleon absorption cross section of cr a b s = 3.6 mb. 



Process 



cr, 



tot 
cxpt. 



n=oo 



n=3 



S+U: 



Pb+Pb: 



Drell-Yan 310 ±10 + 25 nb 449 nb 328 nb 272 nb 

J/tp 7.78 ± 0.04 ± 0.62 fib 12.2 fib 8.38 fjb 6.83 fib 

Drell-Yan 1.49+ 0.01 ± 0.11 fib 1.97 ^b 1.22 fib 0.95 fib 

J /if) 21.9+ 0.02 ± 1.6 fib 44.5 ^b 26.1 fib 19.8 fib 
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